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In view of the recent quest for well-behaved nonlinear ex- 
tensions of the traditional Schrüdinger-von Neumann unitary 
dynamics that could provide fundamental explanations of re- 
cent experimental evidence of loss of quantum coherence at 
the microscopic level, I review the general features and com- 
plete with an additional fundamental ansatz the nonlinear 
quantum (thermo) dynamics I proposed [in my doctoral the- 
sis (MIT, 1981) and the subsequent series of papers cited in 
G.P. Beretta, Found. Phys. 17, 365 (1987)] for a single iso- 
lated indivisible constituent system, such as a single parti- 
cle or spin or atòmic system, or a Bose-Einstein or Fermi- 
Dirac field. The nonlinear dynamics entails a fundamental 
unifying microscopic proof and extension of Onsager's reci- 
procity and Callen's fluctuation-dissipation relations to all 
nonequilibrium states, close and far from thermodynamic 
equilibrium. Together with a self-contained review of the 
main results already proved, I show an explicit geometrical 
construction of the equation of motion from the steepest- 
entropy-ascent ansatz and clarify its exact mathematical and 
conceptual equivalence with the maximal-entropy-production 
variational-principle formulation presented in S. Gheorghiu- 
Svirschevski, Phys. Rev. A 63, 022105 (2001), and I show 
how it can be extended to the case of a composite system 
to obtain the general form of my equation of motion, con- 
sistently with the demanding requirements of strong separa- 
bility. Because both my original geomètric construction and 
its equivalent variational formulation, fix only the direction 
(of steepest entropy ascent) along which the state operator 
is driven by the dissipative term of the equation of motion, 
but leave unspecified the rate at which this happens, I pro- 
pose a possible, well-behaved and intriguing, general closure 
of the dynamics, compatible with the nontrivial requirements 
of strong separability. Based on the time-energy Heisenberg 
uncertainty relation, I derive lower bounds to the internal- 
relaxation-time functionals that determine the rate of entropy 
production. This bounds entail an upper bound to the rate 
of entropy production. Guided by this reasoning, I propose a 
maximal-entropy-production-rate ansatz, by which each indi- 
visible subsystem follows the direction of steepest perceived 
entropy ascent at the highest rate compatible with the time- 
energy uncertainty principle. Subject, of course, to experi- 
mental and further theoretical validation, this closure ansatz 
completes the original nonlinear dynamics with no need of 
new physical constants. 



I. INTRODUCTION 

In my 1981 doctoral thesis 0] and a subsequent se- 
ries of papers (see also Refs. 00), I addressed 
the problem of deriving a well-behaved extension of 
the Schròdinger-von Neumann unitary dynamics meet- 
ing the very demanding set of strict requirements that 
appear to be necessary if one is willing to accept [pd| : 
(1) the "broader quantum kinematics" ansatz that p — 
a unit-trace, nonnegative-definite, hermitian operator on 
the Hilbert space H associated with a system according 
to Standard quantum mechanics — represents a "true" 
quantum state (in the sense analogous to that of the 
wave function of Standard QM) not only when p 2 = p 
but also when p 2 ^ p, even if the system is strictly iso- 
lated and uncorrelated; and (2) that — fc B Tr(plnp) rep- 
resents the "physical" entropy (as opposed to a statisti- 
cal or information-theoretic entropy related to incoherent 
stochastic mixtures of true states). 

The equation of motion I designed and proposed for 
the state operator p of a general system consisting of M 
distinguishablc indivisible subsystems has the form 

^ = -^[H,p}+Ü M (p,H,G l ) , (la) 

M 

(lb) 

d j = lVPj( mn P) J ]±c{ VP ji J , V pj(Hy AVP7Í^) J ]} ' (lc) 

where the notation is defined in Section [x|. 

For a single indivisible system, it reduces to the form 

^ = ~{H, P ]+Í> l (p,H,G i ) , (2a) 

U x (p, H, G t ) = -^[VP D + DÍ Vp\ > (2b) 
D = [«/plnp] xc{V p ItV p H ^ Gi]} , (2c) 

where [^/p ln p]±c{^/pRi} denotes the component of opera- 
tor y/pla. p orthogonal to the linear manifold C{^/pRi} = 
£{^/pI, y/pH, [y/pGi]} spanned by all real linear combi- 
nations of the set of operators {y/pl, y/pH, [y/pGi]} where 
I is the identity operator on the Hilbert space Tí of the 
system, H the usual Hamiltonian operator associated 
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with the system in Standard QM, the GVs are additional 
(not always necessary) hermitian operators commuting 
with H that we call the non-Hamiltonian generators of 
the motion and that depcnd on the structure of the sys- 
tem, such as the z-th-type particle-number operators Ni 
for a Bose-Einstcin or Fermi-Dirac field (in which case 
Ti is a Fock space), or the total momentum component 
operators P, of a free particle (in which case Gheorghiu- 
Svirschevski [Q proved Galilei invariance). 

As I realized in Ref. ||, the internal-relaxation times 
Tjr(p) need not be constants (as initially assumed in Ref. 
P|) but can be any positive functionals of p. However, 
here I note that the condition of separability imposes 



precise restrictions on these functionals. In Section Vil 



I derive a lower bovmd for Tj(p) which corresponds to 
the highest entropy production rate compatible with the 
time-energy Heisenberg uncertainty principle. Taking 
the internal-relaxation times Tj (p) exactly equal to these 
nontrivial lower bounds turns out to be compatible with 
the separability condition, and completes the nonlincar 
dynamics in a general way with no need to assume the 
existence of new physical constants. 

Equation (Q) satisfies all the necessary requirements 
(listed in Appendix |A|) for a self-consistent and well- 
defined extension of the Schròdinger-von Neumann equa- 
tion of motion compatible with thermodynamics require- 
ments, and turns out to have very intriguing mathemat- 
ical properties and consistent physical consequences. 

In this paper, in addition to providing a self-contained 
review of results already proved in previous work, I give a 
more explicit account of the steepest-entropy-ascent geo- 
mètric construction on my nonlinear equation of motion 
for quantum thermodynamics than available in the orig- 
inal papers. As a result, I show that the equivalent vari- 
ational formulation in Ref. |ïl| can be extended to the 
composite system case. I also show that, in order for the 
general form of my equation to satisfy strong separability 
conditions conceptually equivalent to those recently dis- 
cussed by Czachor |15| , the relaxation-time functionals 
must satisfy nontrivial conditions. 

I emphasize that the nonlinearity of the dynamics is 
such that one should not expect that the equation vàlid 
for a singlc clcmcntary (indivisible) constituent of matter 
be vàlid also for a system with an infernal structure, be- 
cause the equation of motion must reflect such structure 
not only through the Hamiltonian operator. However, of 
course, the superoperator Üjf must reduce to Di [see 
Eq. (|94|) below] when the system has only one indivisible 
subsystem. 

For a composite system, my original general formula- 
tion To make the paper self-contained, I briefly review 
the fundamental unifying extensions to all nonequilib- 
rium (dissipative) states of the Onsager reciprocity re- 
lations |ïl| and Callen's fluctuation-dissipation theorem 
JÏ4j that derive from Eq. (|2|) [|| as well as, in general, 
fromEq. (@) §. 

In Appendixes |A] and |b| I discuss a set of criteria and 
the definition of stability of equilibrium that a fundamen- 



tal dynamics must satisfy in order to be compatible with 
the second law. In the other Appendixes, I discuss some 
extensions of the mathematics of Eqs. |ï]) and (0). 



II. THE AUGMENTED STATE DOMAIN ANSATZ 

The fundamental ansatz that the postulates of quan- 
tum mechanics can be successfully supplemented by the 
first and second principies of thermodynamics by assum- 
ing a broader state domain that includes not only p 2 = p 
but also p 2 ^ p state operators, provided that the func- 
tional — fc B Tr(plnp) is taken for the physical entropy, was 
first proposed (without a dynamical law) by Hatsopoulos 
and Gyftopoulos [ïl}| . 

Twenty five years ago, the hypothesis of a state domain 
augmented with respect to that of traditional QM was 
perceived as countercurrent to the prevailing approaches 
to dissipative quantum dynamics within the frameworks 
of statistical, stochastic, phenomenological, information- 
theoretic, chaotic-behavior and bifurcation theories. For 
this reason, the broader quantum kinematics ansatz |IÏ| ] 
and, five years later, my equation of motion have been 
considered initially as unphysical and substantially ig- 
nored, mainly because their motivation appeared to be 
derived from theoretical reasoning only. 

In search for direct experimental evidence, I computed 
the effects of the irreversible atòmic relaxation implicd 
by the nonlinear equation of motion onto some bàsic 
quantum-electrodynamics results on absorption, stimu- 
lated emission, and resonance fluorescence from a single 
two-level atom The results were obtained in the 

near-equilibrium linear limit and, of course, in terms of 
the yet undetermined internal-relaxation-time functional 
r(p) that is part of the equation of motion. To my knowl- 
edge no one has yet attempted to verify these results ex- 
perimentally and estimate r(p). 

The recent new experimental evidence of loss of quan- 
tum coherence jÏ6 17 and the impressive effort dc- 
voted to study nonlinear modifications of the Stan- 
dard Schròdinger equation in the last ten years ||Ï8| , 
finally seem to make more acceptable, if not require, 
P 

Gyftopoulos ansatz. 



the p 2 ^ p augmented state domain Hatsopoulos 



Once the p ^ p ansatz is accepted, the nonlinear 
equation of motion I proposed completes the dynam- 
ics and holds the promise to provide a microscopic-level 
explanation of the recent experimental evidence of loss 
of quantum coherence. It is with this motivation that 
Gheorghiu-Svirschevski Jï^] has "rediscovered" the sim- 
plest form of my equation together with many of its 
known features. Ref. jÏ2| contributes to confirm the va- 
lidity of my equation, including existence and unique- 
ness of solutions, and derives useful results in the near- 
equilibrium linear limit. However, the question of defin- 
ing the form of the relaxation-time functional is left un- 
resolved. 
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III. THE STEEPEST-ENTROPY-ASCENT 
ANSATZ FOR AN INDIVISIBLE SYSTEM 

In Refs. §,§,0, 1 emphasized that Eqs. (0) and (§) have 
an important geomètric interpretation. The Hamiltonian 
tcrm, —í{H,p]/fï, and the dissipative term, T)(p, H,Gi), 
compete with each other in the sense that the first drives 
the state operator toward a unitary motion tangent to 
the local constant entropy surface, whereas the second 
drives it toward the local direction of steepest entropy 
ascent along the surface with constant mean vàlues of 
the generators of the motion. In this section, I show 
more explicitly than in the original papers how Eq. (|^) 
can be constructed directly from the steepest-entropy- 
ascent ansatz. 

Let Tí (dimH < oo) be the Hilbert space and H the 
Hamiltonian operator that are associated with the given 
(indivisible) elementary constituent system in Standard 
QM. For simplicity, we first consider a system composed 
of a single elementary constituent. The generalization to 
AI constituents is given in Section [xj. 

We assume that the (true) quantum states are one- 
to-one with the linear hermitian operators p on 7i with 
Tr(p) = 1 and p > 0. As done in Ref. |2|, we intro- 
duce the square-root state operator y/p obtained from p 
by substituting its eigenvalues with their positive square 
roots fTJ. 

As a first step to force positivity and hermiticity of the 
state operator p we assume that the equation of motion 
may be written as 



dp 
dt 



(3) 



where the operator E, in general non-hermitian, is de- 
fined later. 

We consider the space C{7í) of linear operators on H 
equipped with the real sealar product 



(F\G) = ^Tv(F^G + G^F) 



(4) 



so that for any (time-independent) hermitian R in C(7í) 
the corresponding mean-value state functional, its local 
gradient operator with respect to yfp, and its rate of 
change are 



r(p) = Tr(pR) = (^p\^R) , 
V^r(p) = ^ = v ^R + R y /p, 



(5) 
(6) 

(7) 



Moreover, the entropy state functional, its gradient op- 
erator with respect to y^p, and its rate of change are 

s(p) = -fc B Tr(plnp) = -k B (y/p\^/plií p) , (8) 
ds(p) 



V^ps(p) 



(9) 



dajp) 
dt 



= E 



dp 
dt 
ds{p) \ 
dVP ) 



Tr (^7) + Tr(^lnp) 



E 



dp 
dt 
ds(p) 



(10) 



Now it is easy to see from Eq. (R) that the vàlues of the 
mean functionals ri(p) are time invariant if and only if 
E is orthogonal to yfpRi, for all í, i.e., if it is orthogonal 
to the linear manifold Cl^/pRi} spanned by the set of 
operators {yfpRi} in which we always have yfpRo = ■Jpl 
(to preserve Trp = 1) and <Jp~Ri — ^/pH (to conserve 
energy), plus the additional non-Hamiltonian generators 
of the motion as already discussed. 

It is noteworthy that 



E 



H 



<P)I\ 



(11) 



with c{p) any real functional of p, yields [through Eq. 
(0)] the Hamiltonian part of the equation of motion, and 
is orthogonal to C{yípRi} and to the entropy gradient 
ds(p)/dy / p as well. 

Introducing the notation |^| 

AF = F - Tr(pF)/ , (12) 
(AFAG) = (AGAF) = {^AF^AG) 



= -Tr(p{AF,AG}) , 



(13) 



for F and G hermitian in C(H) and {•, •} the usual an- 
ticommutator, and defining the shortest characteristic 
time associated with the Hamiltonian part of the equa- 
tion of motion by the relation (see Appendix ^) 



^(E H \E H ) , 



(14) 



th{p) 2 

we find, from Eq. ( |ÏT| ) and (y/p\yfp) = 1, 

T H (p) 2 (AHAH) = h 2 /4 - [c(p) + Ti{pH)] 2 T H {pf , 

(15) 

so that only the choice c = — e(p) = —Tr(pH) is compat- 
ible with the Heisenberg time-energy uncertainty rela- 
tion th(p) 2 (AHAH) > h 2 /4 which is then satisfied with 
strict equality. Therefore, 



and 



E H = -^pAH, 



t h (p) 2 (AHAH) = h 2 /A 



(16) 



(17) 



We may note that the foregoing discussion on the 
choice of c(p) is somewhat artificial because any c(p) can- 
cels out in Eq. (||). However, the line of thought becomes 
important when we apply it later to the dissipative part 
of the dynamics. 
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For the general rate of change of ^/p, we let 



E = E + E 



(18) 



and assume E D orthogonal to £{yfp~Ri} but in the direc- 
tion of the entropy gradient operator. We cannot take 
E D dircctly proportional to ds(p)/dyfp as such, because 
the entropy gradient in general has a component along 
C{tJp~Rí\ which would not preserve the mean vàlues of 
the generators of the motion. We must take E D pro- 
portional to the component of the entropy gradient or- 
thogonal to £{^/pRi} 7 namely, for the indivisible system, 



E 



1 



D 4fc B r(p) 
1 



ds( P ) 



_L£{v^fl s } 



(19a) 
(19b) 



d^P 

{ ds{p) 
Ak B r(p) \ d^p 

= ~^(p){^P ln P-^' lnp '\c{ V pR l }) ( 19c ) 

(19d) 



C{^pR l } / 



2t( P ) 
1 



D 



where [ds(p)/dy/p] C { y /-pR i } denotes the projection of 
ds(p)/dy/p onto £{y/pRi} and 



D = Vplnp- I\/Plnp] £{v ^ i} 



(20a) 
(20b) 



y/plnp 


yfpRo ■ ■ ■ yfpRi 




(Vp 


InpIVP-Ro) 


(VP-Ro \^pRa) ■■■ {^pR % \^pRa) •■■ 


(Vp 




(^pR„\^pRi) ■■■ (^pRi\^pRi) ■■■ 




JpAS 


^pAR 1 ••• ^pAR t ••• 


(20c 




{ASAR 1 ) 


(AüiAüi) • • • {AR i AR 1 ) ■ ■ ■ 
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(ASARí) 


(ARíARi) ■ ■ ■ {ARíARí) ■ ■ ■ 













(20d) 

where, for convcnicncc, we dcfincd the entropy operator 

S=-k B Bliip (21a) 

with B obtaincd from p by substituting its nonzero eigen- 
values with unity, i.c., more formally, B = B 2 = PR an/ 3 = 

I — Plierp = P±Kerp, SO that 



B = B 2 , [B,p] = 0, Bp = p. 
The Gram determinant at the denominator, 

{^pRo\Vp~Ro) ■■■ {Vp r í\Vp r o) ■■ 

(ypiïolv^-R.) ••• (^pRi\^pRi) ■■ 



det[{(Aü í> oAü j>0 )}] 

(AiïiAfli) • • • (ARíARí) ■ ■ ■ 
(ARíARí) ■ ■ ■ (AR t ARi) ■ ■ ■ 



(21b) 

(22a) 
(22b) 
(22c) 
(22d) 



is always strictly positive by virtue of the linear indepen- 
dence of the operators in the set {y/pRi}. 

A further, compact expression can be written if we 
choose a set of operators {y/p Ai} which like the set 
{yfpRi} spans L{^fpl, y/pH, [y/pGi]} but, in addition, 
forms an orthonormal set, e.g., obtained from y/pl, y/pH, 
y/pGi by a Gram-Schmidt orthogonalization procedure 
followed by normalization. For example, 



As = 



Ai=I , 
A = AH 

2 y/{AHAH) ' 

{AH AH) Ad- {AH Ad) AH 

\/(A_fíA/í)[{AHAír)(AGiAGi)-(AírAGi) 2 ] 



Then, 



(VpAIVpA) = % . 

D = y/plnp- T,í=i(Vp 1ií p\Vp A í)Vp A í 

r(ÍVpA}) = i • 



(23a) 
(23b) 

(23c) 



(24a) 
(24b) 
(24c) 



It is noteworthy that operators Ai are nonlinear functions 
of p and, therefore, in general, vary with p as it evolves 
with time. 

We finally obtain the following equivalent expressions 
for the rate of entropy change 



1 



ds(p) _ 

dí 4fc B r(p) 



ds(p) 



d^P 



-LC{^pRi} 



ds(p) 
d^P 



4fc B r(p) (E\E) 
4fc B r(p) (E D \E D ] 



k B r(p) 



(ASAS) -^((ASAA)) 2 



r{ P ) 



(D\D) 



_L£{Vpfli}/ 

(25a) 
(25b) 
(25c) 

(25d) 
(25e) 
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fc B r( v ^lnp,{ v ^i?J) 

" t(p) ra^pRi}) 

where the Gram determinant 

(Vp s \Vp r o) (vp^oIvp^o) 



(25f) 
(25g) 



(ASAS) (AfliAS) 
(ASAfli) (AfíiAüi) 

(ASARi) (AR 1 AR i ) 



■ (ARiAS) 

■ (AfiiAfíi) 



(26) 

is also strictly positive except when p satisfies Eq. (|2Ï]). 



IV. SCHRODINGER-VON NEUMANN LÍMIT 

Equation (|2|) reduces to the Schròdinger-von Neumann 
equation of motion, ihp = [H,p], when p 2 = p (thus en- 
tailing the usual unitary hamiltonian dynamics of Stan- 
dard QM), and also when and only when 



Bexp(C) 
P ~ Tr[Bexp(C)] ' 
B 2 = B , [B, C] = , 

or, equivalently, 

Vplnp = J-pC - ^lnTr[Bexp(C)] 

where 

C = -0H + víGí , 



(27a) 
(27b) 

(27c) 

(27d) 



in which case we say that the state is nondissipative, and 
the solution of the equation of motion is 

p{t) = B(t) exp(C)/Tr[S(í) exp(C*)] , (28a) 
B(t) = Ui^BiOp-^t) , U(t) = &qp(-itH/h) , (28b) 

which includes the usual p 2 = p Schròdinger dynamics 
when C is the null operator [[3 — Ui — in Eq. (p7|)] and, 
therefore, Tr(B) = 1. 



When [B(0),H] = the nondissipative state is also 
an equilibrium state [or stationary state, if H = H(t)]. 
When \B(0),H] ^ 0, instead, the unitary evolution in 
Eq. ( p8| ) can be regarded as a (constant entropy) limit 
cycle (or "ridge" as termed in Ref. (Ï2|), which coincides 
with the usual periòdic solutions of the Schròdinger equa- 
tion when B is a one-dimensional projector, B — 
T'í/'': = I]- 

Indeed, Equation (0) preserves the rank and nullity 
of p [|| or, as said in Ref. |p^| , preserves the cardinal- 
ity of the set of nonzero eigenvalues of p. Because the 
equation of motion attracts the state operator towards 
the highest entropy state or limit cycle compatible with 
the initial mean vàlues of the generators of the motion 
and the number of zero eigenvalues, a minor perturbation 
of the state that changes an initially zero eigenvalue to 
an arbitrarily small nonzero value would cause an irre- 
versible departurc of the state towards a different (highcr 
entropy) equilibrium state or limit cycle. Hence, as long 
as there are zero eigenvalues of p, i.e., unlcss B = I, all 
equilibrium states and limit cycles are unstable. With- 
out this conclusion, we could not claim that the equation 
of motion entails the second law of thermodynamics (see 
Appendixes |A| and|B| and the original papers, for further 
discussion of this important point). 

The Hamiltonian operator may be set to vary with timc 
in order to model (adiabatic) energy exchange with the 
surroundings. Thcn the unitary part of the equation of 
motion, —i[H(t), p]/h, may cause energy exchange, but 
not entanglement nor entropy exchange with the sur- 
roundings. The dissipative term, instead, describes in- 
fernal relaxation only and does not contribute to energy 
change in any case. 



V. HIGHEST-ENTROPY EQUILIBRIUM STATES 

The only globally stable equilibrium states of the dy- 
namics generated by Equation || are 



P : 



exp(C) 



Trexp(C) Trexp^iT + Vi Gi) 



(29) 



The definition of global stability, which is stronger than 
Lyapunov or local stability and is required by the second 
law l20U2Ï|] , is given and discussed in Appendix The 



ffland relied on a technical 
201 and later found to be 



proof was discussed in Refs. (| 
conjecture discussed in Ref. 
proved by Hiai, Ohya, and Tsukada p2[ . 

For a system with generators of the motion /, H, Ni, 
Eq. (|2^) represents in general the grand canonical ther- 
modynamic equilibrium states. It reduces to the canon- 
ical equilibrium states if = n^J for all i's (to the mi- 
crocanonical if also H = el) . 

As is well known, states given by Eq. (|2^) are solutions 
of the constrained maximization problem 



max s(p) 

subject to »"j(p) = Ti and p > , 



(30a) 
(30b) 



5 



where s(p) = -fc B Tr(plnp), n(p) = Tr(pJ), r 2 (p) = 
Tr(pií), ri(p) = Tr(pGi), n = 1, r 2 and rj's given. 

The inequality constraint p > can be eliminated re- 
casting the problem in terms of ^/~p, 

max s(y/p) (31a) 
subject to n^/p) = r t , (31b) 

where s(^fp) = -fc B Tr[( v /p) 2 lní^) 2 ] = s(p) and 
ri(y/p) = Tr[(y/p) 2 Ri] — ri(p). The method of Lagrange 
multipliers then gives the condition (necessary but not 
sufhcient) 



ds(p) 



5> 



dn{p) 



= 



(32) 



9 Vp i d \fp 

which, using Eqs. (||) and ([}]), becomes 

- 2k By /phip-2k BV fp-Y, X i (VpRí + Ri \ÍP) =° • 

i 

(33) 

It is noteworthy that Eq. ( |33| ) is satisfie d wi th obvi- 
ous identification of multipliers by Eq. (27c), where 
Ei ^iyfpRi commutes with ^fphip, i.e., with p. Thcrc- 
fore, each nondissipative state satisfies the necessary con- 
dition (|3^) although it is not a solution of the maximiza- 
tion problem (J3l|) unless B = I. 



VI. ONSAGER'S RECIPROCITY AND 
CALLEN'S FLUCTUATION-DISSIPATION 
NONEQUILIBRIUM RELATIONS 

First, I introduce a general representation of state op- 
erators particularly useful for representing nonequilib- 
rium states. Any p can be written as [gj 



P = 



BeM-EjfjXj) 
TrBexp(-V/ 3 X 3 ) 



(34) 



where the hermitian operators in the set {/, Xj} span 
the real space Ch{Ti) of linear hermitian operators on Tl, 
and B is the usual idempotent operator defined in Eq. 
(211). Indeed, the operator B\np is always well-defined 



and belongs to Ch(Ti), so that 

B Inp = -f I - fi X i = -foBI - Ç fjBXj (35) 

3 3 

where the second equality follows from B 2 = B. Thcrc- 
fore, 

yfp lnp = -foy/p - Ej fi \ÍP X 3 . (36) 

x J (p)=Tr(pX 3 ) , (37) 

s(p) = fc B /o + Ej /i (p) i (38) 

D = -Ej fi [Vp x 3 - EL·i(Vp x 3\Vp a ^)VpM (39) 



where, for simplicity, we use the expression (24b) for D 
in terms of the orthonomal set {y/p Ai} [an equivalent 
expression in terms of Gram determinants is readily ob- 
tained using Eq. Q i n Eq. (|(]c)]. 

Notice from Eq. (p8h that the partial derivative 



kafj — 



ds{p) 



dxj (p) 



(40) 



,(/->) 



may be interpreted as the generalized afhnity or force 
representing the entropy change that corresponds to an 
independent change in the mean value of the linear ob- 



servable Xj. 



Next, we define the dissipative rate of change of the 
linear mean- value functional associated with operator Xj, 



Dí 



Trp^p, H, Gí)Xí] =2(E D \ JpX t ) . (41) 



Using E D = -D/2r{p) [Eq. pd| )] and Eq. @ yiclds 
the linear interrelations between dissipative rates and 
generalized afhnities, 



Dí 



^2fi L a(p) 



(42) 



where the coefhcients may be interpreted as generalized 
(Onsager) dissipative conductivities and are nonlincar 
functions of p, 

L l0 {p) = L^p) = -L-K^pX^X,) 

~ Y,Li(VP~ X i\VP A k)(VP A k\VP X i)] W 
[(AXíAXj) -Efc=i (AI,·AA) t )(AA)tAI' J ·)] (43b) 

([\/P X i]±C{^pR i }\[</P X i]±C{^R i }) (43c) 



-00 



r{p) 





(AXiAXj) 


(ARxAXj) ■ ■ ■ (ARkAXj) ■ ■ ■ 






(AfliAiïi) • • • (ARkARi) ■ ■ ■ 






(AR 1 AR k ) ■ ■ ■ (AR k AR k ) ■ ■ ■ 
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r(p) 




mvpRk}) 



(43d) 



This equations are at the same time a proof and a gener- 
alization of Onsager's reciprocity relations ]Ï3| extended 
to all nonequilibrium (dissipative) states. Moreover, they 
give explícit expressions for the nonlincar dependence of 
the dissipative conductivities Lij(p) on the state oper- 
ator p, the generators of the motion and the internal- 
relaxation-time r(p). 

Using Eqs. (]2^) and ([36|), the rate of entropy change 
may be rewritten as a quadratic form of the generalized 
afhnities, 



G 



ds(p) 
àt 



(44) 



The symmetric matrix [{L !;| -(p)}] is a Gram matrix [see 
Eq. (43c)] and, as such, it is nonegative definite, i.e., its 



determinant 

det[{Lij(p)}] > , 
with strict positivity only if all operators 



(45) 



(46) 



are linearly independent, in which case Eq. ( [42| ) may be 
solved to yield 



~Dxi(p) 
Dt 



(47) 



and the rate of entropy change can be written as a 
quadratic form of the dissipative rates 



Mp) _ / CB ^^2r.i(p)- DXl ^ D:Ej ^' ) 



dí 



Dí 



Dt 



(48) 



Equations (43) are also a proof and generalization of 
Callen's fluctuation-dissipation theorem extended to 
all nonequilibrium (dissipative) states. Indeed, we intèr- 
pret (AXiAXj) as the codispersion (covariance) of mea- 
surement results of observables Xj and Xj when the sys- 
tem is in state p and (AXíAXí) as the dispersion (or fiuc- 
tuations) of measurement results of observable Xí. As in 
Callen's fluctuation-dissipation theorem, the expressions 
in Eq. ( pÏ3| ) relate codispersions with generalized conduc- 
tivities. 

A judicious choice of the set {^/pXj} may greatly sim- 
plify these relations. In particular, if {y/pXj} is any or- 
thogonal extension of the orthonormal subset {yfpAi}, in 
the sense that {^fpXj\^fpAi) = for all Aj's and j > a, 
then for all i and j > a, 



Lij(p) 



-00 



(AXiAXj 



(49) 



which relates directly [through the internal-relaxation 
time t(p)\ covariances and fluctuations in the observables 
Xi with their associated dissipative conductivities. 

Onsager's 13 and Callen's [Q theorems are keystones 
of our understanding of irreversibility. Indeed, the proof 
I proposed emerges not from the analysis of the Hamilto- 
nian term of the equation of motion supplemented with 
the so-called assumption of microscopic reversibility, but 
from the dissipative term of the equation of motion, i.e., 
the only term responsible for irreversibility, with no ad- 
ditional assumptions. 

Onsager's result JÏ3[ was obtained from empirical ob- 
servations on nonequilibrium phenomena very close to 
stable thermodynamic equilibrium, so that the list of X^'s 
was indeed very short, and the result vàlid only for a 



limited class of states. Our result generalizes the validity 
of Onsager's reciprocity relations to all nonequilibrium 
states, close and far from stable thermodynamic equilib- 
rium. Of course, the price we have to pay to describe 
nonequilibrium states far from stable equilibrium is that 
we must use a much larger, possibly infinite list of Ai's 
in Eq. m. 



VII. INTERNAL RELAXATION TIME VS 
TIME— ENERGY HEISENBERG UNCERTAINTY: 
MAXIMAL-ENTROPY-PRODUCTION-RATE 
ANSATZ 

By analogy with what done for the Hamiltonian char- 
acteristic time th(p), we define the shortest characteris- 
tic time of the dissipative rate of change of y/p by the 
relation (see Appendix |(]) 



1 



td{pY 



From Eqs. (|Ï9| ) and (gC 
from Eq. (Eq)] , we find 

td( P ) 2 = klr( P f 



±(E D \E D ) (50) 
i and (y/p\y/p) = 1 [or, directly, 



r(vpS,{VM) (^1^) 



(51) 



From this relation we may extract a possible ansatz 
on the functional form or at least a lower bound for the 
internal-relaxation-time functional r(p) in Eq. (|^) by as- 
suming that also td(p), like th(p), should satisfy the 
Heisenberg uncertainty relation 



T D (p) 2 (AHAH) > h 2 /4 



This implics 



-(P) 2 > 



h 2 r(^pS,{^pR t }) 



h 2 (D\D) 



4k 2 (AH AH) T({^pR t }) A(AHAH) 
or, using Eqs. (||) and (fy|), 

h 2 

T{P) - A(AHAH) ^^ fifjLij{p) 



(52) 



(53) 



(54) 



Relation ( |53| ) implies an upper bound to the rate of en- 
tropy production, 



ds(p) 2 (AHAH)T(JpS : {VpRi}) . . 

By analogy with the Hamiltonian time, a possible ansatz 
that I propose should be carefully considered is that also 
td(p) and, therefore, r(p) be equal exactly to the lower 
bound [i.e., strict equality in Eqs. (fÏ3|), ( p^ ) and (p5[)], 
corresponding to the maximal entropy production rate 
compatible with the time-energy uncertainty relation. 
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This corresponds to a truly maximal entropy production 
dynamics. The square-root state operator ^fp not only 
moves in the direction of steepest entropy ascent, but it 
does so at the highest rate compatible with the uncer- 
tainty principle. 

This ansatz is very intriguing, conceptually appeal- 
ing and fraught with far-reaching physical conscquences. 
But, like any other hypotheses on the functional forní 
of r(p), it should be verified against known behavior in 
worked-out specific cases, such as the examples consid- 
ered in Refs. ||,||,[Ï2| and the variety of physical prob- 
lcms for which nonlinear modifications of the Schrodinger 
equation have been deemed necessary but have not yet 
been resolved [ ^8| . 

If proved vàlid, this ansatz would complete the dy- 
namics without need of additional physical constants and 
would provide a fundamental microscopic foundation of 
the macroscopic observation that Nature always evolves 
at the fastest possible rate along the most direct path 
towards maximum entropy compatible with the system's 
structure and external constraints. Interpreted in this 
way, my nonlinear dynamics would provide a unifying 
fundamental microscopic foundation of all phenomeno- 
logical theories of irreversible processes advanced in the 
last seventy years after the pioneering work of Onsager. 



VIII. VARIATIONAL PRINCIPLE 
FORMULATION 

Following Gheorghiu-Svirschevski |ï^| , the direction of 
steepest-entropy-ascent can also be found by considering 
the constrained maximization problem 



maX d!^ 
subject to -^-(E D ) 



and (E D \E D ) = 



(56a) 
(56b) 



where fi and r|, are given real functionals of p. The nec- 
essary solving condition in terms of Lagrange multipliers 
is 



d ds 
dE D dt 



-5> 



d àrj 
dE D dí 



d 



dE T 



(E D \E D ) = 0. (57) 



Using Eqs. (0) and (JÏC|) , it becomes 



(58) 



where the multipliers must be determined by substitution 
in the system of constraint equations. 

It easy to verify that our expression of E D — —D/2t(p) 
m Eqs. (^ and © is an explicit solution of Eq. ( |58| ) in 
the case fi = and T£> given by the explicit expression 
in Eq. (pT"l). See Appendix ^ for f t ^ 0. 



IX. BLOCH NEAR-EQUILIBRIUM LINEAR 
LÍMIT 

In Ref. Q we have shown that in the neighbourhood 
of each equilibrium state p c {given by Eq. (|27]) with the 
additional condition [B, H] = 0}, if we assume r(p) con- 
stant (e.g., r = r(p c ) ^ 0), the state operators in the 
subsct with B = B e and Tr(pRi) = Tr(p c Ri) (i.e., with 
the same nullity, rank and mean vàlues of the generators 
as the equilibrium state) obey a linearized form of the 
equation of motion that has the form of a Bloch relax- 
ation equation, 



dp 
dí 



■[H,p] 



P~ Pc 



so that the solution of the equation of motion is 
p(t) -» expi-t/r^U^piOp-'it) 

+ [1 - exp(-í/T e )] p e (0) , 



(59) 



(60) 



with U(t) given by Eq. ( 28b ) {note that if there are no 
non-Hamiltonian generators the condition [B, C] =0 im- 
plies [p,H] = and, therefore, U^p^U^ 1 ^) = p(0)}. 

Gheorghiu-Svirschevski p2| , derived a linearized form 
for the case with no non-Hamiltonian generators vàlid 
also when [p, H] ^ (but assuming r — > r c ^ 0), and 
showed that it yields a generalized Fokker-Planck equa- 
tion for a free particle, and a Langevin equation for a 
harmònic oscillator. 

It is however noteworthy that such linearized limit be- 
havior should be revised if we take r(p) equal to the 
lower bound imposed by the time-energy uncertainty re- 
lation [Eq. (53) with strict equality], because then, if 
(AH AH) yí 0, r(p) -> as p -> p c . 



X. DIVISIBLE SYSTEMS DYNAMICS. 
DEFINITIONS AND NOTATION 

The composition of the system is embedded in the 
structure of the Hilbert space as a direct product of the 
subspaces associated with the individual elementary con- 
stituent subsystems, as well as in the form of the Hamil- 
tonian operator. 

In this section, we consider a system composed of 
M distinguishable and indivisible elementary constituent 
subsystems. For example, each subsystem may be a 
different elementary particle or a Fermi-Dirac or Bose- 
Einstein field (in which case the corresponding TÍ J is a 
Fock space). The Hilbert space is 



H = H^H 2 ® ■ ■ ■ 
and the Hamiltonian operator 

AI 



M 



h = J2 H .j®h + v . 



(61) 



(62) 



j=i 
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where Hj is the Hamiltonian operator on H J associated 
with the J-th subsystem when isolated and V (on 7Ï) the 
interaction Hamiltonian among the M subsystems. 

For convenience, we denote by Ji J the direct product 
of the Hilbert spaces of all subsystems except the J-th 
one, so that the Hilbert space of the overall system 

H = H J ®H 7 (63) 

and the identity operator / = F^I-j. 

The subdivision into elementary constituents, consid- 
ered as indivisible, and reflected by the structure of the 
Hilbert space H. as a direct product of subspaces, is par- 
ticularly important because it defines the level of descrip- 
tion of the system and specifies its elementary structure. 
This determines also the structure of the nonlinear dy- 
namical law I proposed. In a sense, this is a price we have 
to pay in order to free ourselves from the assumption of 
lincarity of the law of motion. 

In other words, the form of the superoperators T>m and 
E>i arc different dcpcnding on whcthcr the system is or is 
not subdivisible into indivisible subsystems. For exam- 
ple, consistently with the conditions listed in Appendix 
, we request that the superoperator T)m satisfies the 
strong separability conditions |p|,[Ï5| 



D M {pa®Pb , H A ®I B +I A ®H B , G iA ®I B + Ia®G íb ) 

= &m a {pa,Ha,Gía)®Pb+Pa®& Mb (p b ,H b ,G íb ), (64a) 

Tt b [Í)m(p, Ha®Ib+Ia®H b ,Gía®Ib+Ia®G íb )} 

= F MA (p,H A ,G lA ) for any p, (64b) 

T± A \P M {p,H A ®I B +I A ®H B ,G iA ®I B +I A ®G iB j\ 

= F Mb {p,H b ,G íb ) for any p, (64c) 

where subsystems A and B are obtained by parti- 
tioning the set of constituents 1, 2, . . . , M into two 
disjoint subsets of M A and M B constituents, respec- 
tively {Ma + M B = M). Note that, of course, if 
Pa®Pb then F Ma {p : H a ,G íA ) = t> MA {p A ,H A ,G lA ) and 

Fm b {p,H b ,Gí B ) = T>m b (pb,Hb,Gíb)- 

It is noteworthy that trying to apply the same con- 
ditions to superoperator Di would be physically mean- 
ingless, because if the system were subdivisible then Di 
would not be i ts di ssipati ve e volution superoperator. 

Conditions ( |64rj| ) and (34c) prevent locality problems 
because they guarantee that changes of the Hamiltonian 
(or the other generators of the motion) in one of two non- 
interacting subsystems cannot aífect the mean vàlues of 
local observables of the other subsystem. For example, 
assume that subsystems A and B are correlated but not 
interacting. We may switch on a measurement appara- 
tus within B and therefore alter the Hamiltonian H B . 
By virtue of Eqs. ( |64bj ) and (34c), the rate of change 
of the reduced state operator p A — Tr B (p) does not de- 
pend on H B and, therefore, all functionals of p A (local 
observables) remain unaffected by the change in B, i.e., 
no faster-than-light communication can oceur between B 
and A (of course, if we exclude the projection postulate). 



In addition, we must consider the following additional 
non-trivial conditions of separate energy conservation of 
noninteracting subsystems for any p (see Appendix [A|, 
Condition 6) 

Tr[(H A ®I B ) t> M {p, H A <Z>I B +I A ®H B , G t )] = , (65a) 
Tr[{I A ®H B ) t> M {p, H A (Z>I B +I A ®H B , G t )] = , (65b) 

and of separate entropy nondecrease for uncorrelated 
(possibly interacting) subsystems (see Appendix |A|, Con- 
dition 7) 

Tr[(S A ®I B ) Í) m {pa®Pb,H, d)] > , (66a) 
Tr[(I A ®S B ) I) M {pA®PB,H, Gi)} > , (66b) 

where Sa — ~k B B A hip A , S B = —k B B B ln p B . 

For each type of particle in the system, we may write 
without loss of generality the number-of-particles-of-ï-th- 
type operator associated with the system as p3| 



N > = Y. N l· J® I J for i = 1,2, 



(67) 



where N iJ denotes the number-of-particles-of-i-th-type 
operator associated with the J-th subsystem. 

In general, we assume that the set of linear hermitian 
operators /, H , Gi on 7i, always including / and H, are 
the generators of the motion of the composite system. 
For example, the list of non-Hamiltonian generators Gi 
may coincide with that of the number operators defined 
in Eq. ©. 

For each subsystem J we denote by £(H J ) the space 

of linear operators on 7i J equipped with the real sealar 
product 



(FtjGj + CfijFj 



(68) 



For a given state operator p on 7í, given linear hermi- 
tian F and G on TL, and each subsystem J, in addition to 



that already defined in Eqs. (^), (|Ï2 
the following convenient notation |2 



and (O), we define 



P j - Tr 7 (p) , (69) 

Pj = Tr. ; (p), (70) 

(F) J = TrjKl^pjjF] , (71) 
(AFAG) 3 = {^ p - J {AF) J \^p- J {AG) J ) J 

= lTT J {p J {{AFy,{AG) J }), (72) 



Tr j over 



where Trj denotes the partial trace over 7i J 

TÍ J , Pj is the reduced state operator of elementary sub- 
system J and pj that of the composite of all other sub- 
systems. 

In view of the special role they play in the equation of 
motion, we intèrpret the operators {H) J and {S) J , 
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(H) J = TvjKIj^H] , (73) 
(S) J = Tr jKIj^S] , (74) 

respectively, as "internal perception" operators repre- 
senting the way the overall-system Hamiltonian and en- 
tropy operators are "felt" locally within the J-th con- 
stituent subsystem. 

For a given pj, and given linear hermitian Fj and G j 

on TÍ J , we further define the notation 



{*F J l·G J ) J = (JpjAF J \JpjAG J ) J 

= ±Tr j(p j{AFj, AG j}) , (75) 

VP] ■ ( 76 ) 

Bj , (77) 

S J = -k B B J ]jip J , (78) 



where the operators y/p~J, Bj and Sj are obtained by sub- 
stituting each nonzero eigenvalue pi of Pj with its positive 

square root, unity and — k B \npi. In general, Sj ^ (S) J . 

The entropy is defined for a subsystem J only if it is 
not correlated with the other subsystems, i.e., if p can be 
written as 

p = pj®pj ; (79) 

then the subsystem entropy is given by the nonlinear 
statc functional of the reduced state operator, 

s j(p.j) = Tr ÀP.j s .j) = -k^jÍPj^Pj) 

= Tr\p(Sj®Ij)] = [y/P^y/PjSjij , (80) 

S = Sj^Ij + Ij^Sj , (81) 
S (P) = s j(Pj) + s j(Pj) > (82) 

and we also have (AS) J — ASj. 

When subsystem J is correlated, instead, its entropy 
is not defined; however, the functional 

s J (p) = Tr j[pj(Sy] = (VP]\VPl( s y)j (83) 

may be interpreted as the subsystem's local "internal per- 
ception" of the overall-system entropy. Only when J is 
uncorrelated, s J (p) = s(p) = s j{pj) + s—(pj). 

The energy is defined for a subsystem J only if it is 
not interacting with the other subsystems, i.e., if H can 
be written as 

H = H,®I- + Ij®Hj ; (84) 
then it is given by the mean-value state functional 

ej(Pj) = ^j(PjHj) = (Vp]\Vp] H j)j > (85) 

and we also have (AH) J = AH j . 

When subsystem J is interacting, instead, its energy 
is not defined, however the functional 



e J (p) = Trilli)'] = (VP}\VPj( H ))j (86) 

may be interpreted as the subsystem's local "internal per- 
ception" of the overall-system energy. Only when J is 
non-interacting, e J (p) = e(p) = e 7 (p 7 ) + ej(p-j). 

The number-of-particles-of-ï-th-type of the overall sys- 
tem in state p and of each subsystem J are given by 
the mean-value state functionals Tii(p) — Tr (pNi) and 

«<» = ^Ip^íj^j)] = n u(Pj) = ^jbjgu) = 
{y/Pj\y/Pj^ij)ji respectively, and from Eq. (p7|) we 

clearly have Tii(p) = y^j— i n iJ (pj). It is noteworthy 
that, for any p, 

(^®^) J = ^- ( 87 ) 

For generality, we assume that the generators of the 
motion on Ti are /, H and G;, with [Gi,H] = but 
not necessarily with the separated structure of the num- 
ber operators. Moreover, for each state p and each con- 
stituent J, we denote by {R u , i = 0, 1, 2, ... , z J } a set 
of hermitian operators in C(Tl) such that the operators 
in the set { y /pj(R iJ ) J }, where 

(Bu) J = T±j[(Ij®Pj) R íj\ > ( 88 ) 

are linearly independent and span the linear manifold 
generated by the operators y/pjl 7 , ^/p~j{H) J , ^J~p~j{Gi) J . 
If the latter operators are linearly independent, for sub- 
system J, then the set {R^} may be chosen to coincide 
with the generators of the motion. If they are not inde- 
pendent, then it could be a smaller subset (in which it is 
convenient, though not necessary, to maintain R QJ = I 
and, if possible, R 1J — H). In any case, we call the oper- 
ators in the set {-R ii7 } the (instantaneous) "generators of 
the motion of subsystem J" . The structure of the dissi- 
pative term Du is invariant under transformation from 
one set {R^} to any other {R' u } with the same defining 
properties. 

For each instantaneous generator R iJ of subsystem J, 
we define the local "internal perception" mean functional 

rf» = Tr jlpjiR^y] = (^1 VPÏ&ijYh ■ ( 89 ) 
We denote by 

£{Vp]( r íj) j } = £{Vp~jIji sfp-j{H) J , Vp7(g,) j } (90) 

the linear span of the operators {^/Pj{Ríj) J } OI , that 
is the same, the linear span of the operators ^fpjlj, 

^p-j{Hy, ^pj(G t y. 

By definition of the set {R^}, the Gram determinant 

r ({Vp]( R t j) J }) = dct [{(Vp7(^j) J lVp7^j) J )j}] 

= det[{(AR u AR jJ ) J }] (91) 
is always strictly positive. 
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The set {R^} can be conveniently chosen so that 
{^/WjiRi j) J } is (instantaneously) an orthonormal set, in 
which case we denote it by { v /~pj(A iJ ) J } and 

r J ({ y /pjíA iJ y}) = i . 



(92) 
(93) 



XI. THE STEEPEST-ENTROPY-ASCENT 
ANSATZ FOR A COMPOSITE SYSTEM 



and (|Ï8|), we now 

(94) 



Maintaining the validity of Eqs. 
assume 

M 
J=l 

as a first step to guarantee the strong separability condi- 
tion (|64| ) [see Appendix [g| for a discussion related to the 
form of Eq. @]. 

The second step is to make sure that E DJ1 which in 
general may be a function of the operators p, H and 
Gi on 7í, reduces to a function of the operators p } , H 7 

and G.j on 7í J only, whenever the constituent is, at the 
same time, uncorrelated (p — Pj ®p-j), non-interacting 
(H + 1 j®H-j) and not coupled through the non- 
Hamiltonian generators {G iJ ®I-j+ 1 j®G rj) . 

To preserve the formal analogy with the notation for 
the indivisible system, we define the local "partial gradi- 
ent" operators 



ds J (p) 



(sy 



~Pj{ R íjY + {Ru) J VP~J > (95) 



(96) 



Now, it is easy to show that the structure of dp/dt 
assumed with Eqs. (|), © and (|J) yields @] 



dt 

dt 



M 



dt 



, (97) 



-fc B Tr(^)-fc B Tr(^51np) 



M 



Eh2fc B (E dj \^Ij)j 



+ 2 ( e dj\Vp](S) j 




(98) 



Finally, we assume that each E DJ is: (A) orthogonal to 
C{^rp~j(R iJ ) J }, so that all the rates of change in Eq. ( p?| ) 
and the first term in the rate of entropy change in Eq. 
( p8| ) are zero; and (B) in the direction of the local par- 
tial gradient of the "internally perceived" overall-system 
entropy functional, s J (p), i.e., we take 



1 



ds J { P ) 
d VPj 



J(5)J 



(99a) 



(99b) 



where, similarly to what done for the indivisible system, 
ÍTjiBlnpY 

j j - L· -WP. J y n i.j, 
fp-AB\iíp\ J 



D .,= 



= ^p](B\iip) J - \^p-j{B\np) J 



(100a) 
(100b) 









■ vpj( AR uy ■■■ 




<ASAiï 17 ) J 




■ (AR.jAR^y ■■■ 




{ASAR i r ) J 


(AR^AR.^ ■ 


■ (AR.jAR^y ■■■ 


1 














J }) 



(100c) 



All the results found for the single constituent extend 
to the composite system in a straightforward way 
The rate of entropy change 



dt 



AI 



/t'n 



AI 



£> b t» (E DJ \E_ 



DJ 1 J 



(101a) 
(101b) 



j=i 

AI 



i r( V p7(5) J ,{yp7(^ J ) J }) 

The dynamics reduces to the Schròdinger-von Neu- 
mann unitary Hamiltonian dynamics when 



fp-j{B\u P y = s f P - J {cy , 
c = \j R u > 



(102a) 
(102b) 



for each J and, therefore, the state is nondissipative 
{equilibrium if [B,H] = 0, limit cycle if [B, H] ^ 0, 
including the case of pure-state Standard QM when 
Tr(S) = l}. 

The maximum entropy thermodynamic equilibrium 
states are given by Eq. (E9T) . 
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Onsager's reciprocity relations follow again from Eq. 
( |35| ) for the operator B ln p. Equations (EJ) and jj^f ) 
are still vàlid [of course, with Dji in Eq. (ftl|)l, with the 
dissipative conductivities given by 



(103a) 



J=i 



= Ljiip) 





(AA' 1 AA'.,) J 




■ (AR^AXj)' ■ ■ ■ 




(AX.Aüi)'' 


(AR 1 AR 1 ) J • 


■ (AR k AR 1 ) J • • • 






(AR!AR k } J ■ 


• (AR k AR k ) J ■ ■ ■ 


1 









r(p) 



n{vp-j{R tJ ) J }) 



(103b) 



Callen's fluctuation-dissi pati on relations, implied by 
the explicit structure [Eq. ( |103| )] of the Onsager dissipa- 
tive conductivities, are greatly simplified if we choose the 
orthonormal set {^fpj(A i } ) J } instead of { y J~fTj{R iJ ) J } 
and the set {^p~j{X.j) J } to be an orthogonal extension 



(104) 



The conjecture about the time-energy Heisenberg un- 
ccrtainty relation in Section VII can be extended as well 



by assuming that each E DJ gives rise to a characteristic 
time T DJ (p) defined by (see Appendix|c|) 



*( E Dj\ E Dj) > 



(105) 



which should independently satisfy the uncertainty rela- 
tion 



r DJ (p) 2 (AH AH) J > tf/A 



(106) 



Therefore, using Eqs. (J99|) and (|100| ), we obtain the lower 
bounds to the intcrnal relaxation times, 



,(P) 2 > 



h 2 r(^py(sy,{^pj(R u y}) _h 2 (D J \D J ) J 



J " ' 4k2 z(AHAH) J T({ V pj(R ij y}) A(AHAH) J 

(107) 



or, using, Eqs. ( pL0l|) and (^4j), 

t j(p) > 



1 AU All , 

and the upper bound to the entropy production rate 



(108) 



ds(p) ^2 



(AHAH) J n V p ] (sy,{ y /p](R u y}) 



mvp-MuY}) 



(109) 



Taking strict equality in each Eq. (107) yields the max- 
imum entropy production rate compatible with the un- 
certainty relations. 

The variational formulation can be found, after assum- 
ing the structure of Eq. (|94|), by solving the constrained 
maximization problem 



ds 

max — {{E DJ }) 



dt 

subject to 



dr 



dt 



(110a) 

(Í e dj}) = r u and (E DJ \E DJ ) = , 



D. 7 

(110b) 

where r iJ and t? are given real functionals of p. By 
virtue of Eqs. ( |97| ) and (p8[), the necessary solving con- 
ditions in terms of Lagrange multipliers for each E D 7 
become 



ds J {p) 



2£V 



P](R 



iJ) 



XtE 



(sy 



(111) 



clearly verified by the E D ys given in Eqs. (|99j) and (100). 

If two subsystems A and B are non-interacting but in 
correlated states, the reduced state operators obey the 
equations 



àpA 
dt 



àpB 
dt 



--t[Ha,Pa] 



1 



M 
JçA 



ypJDj+D\jpj\®(p A }j , (112a) 



1 



M 

JeB 



[VpJ D J+ D \Vp]]®(Pb)j > (H2b) 



where (p A )j = Tt j(p A ) and (pb)j = TtJpb)- 

The strong separability conditions (p4) are satisfied 
provided the internal-relaxation-time functionals are ei- 
ther all constants or satisfy the following set of nontrivial 
conditions. 



(113a) 



Tj(pa®pb,Ha®Ib + Ia®H b , G iA ®Ib+Ia®G iB ) 

Tj(pa,H a ,Gía) for J e A 
Tj(pb,Hb,Gíb) for J e B ' 

Tj(p, H a ®I b +Ia®Hb,G íA ®Ib+Ia®Gíb) 

Tj'{p,H A ,G iA ) for J e A , . 

Tj'(p,H B ,G iB ) for JeB for any p. (113b) 
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It is noteworthy that Conditions 113 are satisfied by the 



maximal-entropy-production-rate ansat z, i.e ., if we as- 
sume that each Tj(p) is given by Eq. (107) with strict 
equality. Inde ed, if A and B arc nonintcracting, the 
structure [Eq. ( |100| c)] of each operator D j with J e A is 
such that any dependence on Hb cancels out, moreover 

(AHAH) J = (AH A AH A ) J Thus, any dependence on 
Hb cancels out in Eq. ( |112| a) an d, si milarly, any depen- 
dence on Ha cancels out in Eq. (112 b). 

It is interesting, however, that dp^/dí in general may 
depend not only on the "local" (reduced) state operator 
Pa but also on the overall state p through the operators 
(B\np) J (with J G A), thus determining a collective- 
behavior effect on the local dynamics originating from 
existing correlations. In fact, operators D j are in terms 

of (AS) J which may differ, when the subsystems are cor- 
related, from operators (ASa) j ■ In other words, the lack 
of interactions between two subsystems does guarantee 
that the energy of each subsystem is conserved [Eqs. (65)] 
and that the reduced dynamics of each subsystem is local, 
in the sense of independent of the Hamiltonian operator 
(and the other generators) of the other subsystem, but 
does not necessarily imply that each subsystem evolves 
independently of existing correlations with the other sub- 
system. 

Regarding the time evolution of correlations, in Ref. ||] 
we defmed the correlation functional between two subsys- 
tems A and B, 



ctab(p) = Tr(plnp) - Tr A (pA lnp^) - Ttb(pb lnps) , 

(114) 

which is nonnnegative definite in general, and zero only 
if p = pa®Pb- The rate of change of the correlation can 
be written as 



à<?AB{p) 
dí 



&AB\H — &AB\D 



(115) 



Based on our understanding of the equation of motion, we 
conjectured that &ab\d should always be nonnegative, if 
it is true that the dissipative term can only destroy cor- 
relations between subsystems, but cannot create them. 
However, this conjecture remains to be proved. The cor- 
responding entropy inequality to my knowledge has not 
yet been studied. 



XII. CONCLUSIONS 

I reviewed most previous results in Refs. [|]-|| on the 
well-behaved nonlinear equation of motion I proposcd 
for quantum thermodynamics that entails the second 
law and nonequilibrium steepest-entropy-ascent dynam- 
ics with Onsager's reciprocity and Callen's fluctuation- 
dissipation relations. Together with the variational prin- 
cipio formulation derived in Ref. |12| and the observation 



in Refs. [ [lipq ] that only the functional — k B Tr(phip) 
can represent the physical entropy, we may conclude 
that all the results so far confirm that my equation 
of motion has all the necessary features to provide a 
sclf-consistent and conceptually-sound resolution of the 
century-old dilcmma on the nature of entropy and ir- 
reversibility, alternative to Boltzmann's statistical ap- 
proach and vàlid also for systems with few degrees of 
freedom. 

The nonlinear dynamics encompasses within a uni- 
fied framework all the successful results of quantum me- 
chanics, cquilibrium and nonequilibrium thermodynam- 
ics. It also holds the promise to provide a fundamen- 
tal framework within which to address the currently un- 
explained evidence on loss of quantum coherence, to 
design new fundamental experiments, to examine new 
applications on the lines of those developed in Refs. 
P,p|JÏ^|, and to further investigate the dependence of 
thc intcrnal-relaxation-time functional on the state op- 
erator and physical constants, as well as possibly verify 
the ansatz proposed in this paper by which each indivisi- 
ble subsystem follows the direction of steepest perceived 
entropy ascent at the highest rate compatible with the 
time-energy uncertainty principle. 

The equation of motion satisfies the set of conditions 
discussed in Appendix^ and, therefore, preserves most of 
the traditional conceptual keystones of physical thought, 
including a strongest form of the non-relativistic principle 
of causality, by which future states of a strictly isolated 
system should unfold deterministically from initial states 
along smooth unique trajectòries in state domain defined 
for all times (future as well as past |pq| ). Interestingly, 
whilc the maximum entropy states are attractors in for- 
ward time, the unitary-solutions boundary limit cyclcs of 
Standard quantum mcchanics are attractors in backward 
time. 

As pointed out by Onsager and Machlup the 
fluctuation-dissipation relations cannot be derived in any 
rigorous way from the traditional Hamiltonian dynamical 
principies, unless these are complemented by some addi- 
tional postulate closely related to the additional princi- 
pies, assumptions, or approximations needed to derive 
the Onsager reciprocity relations. This is sometimes re- 
ferred to as the irreversibility paradox. In other words, in 
order to infer any feature of irreversibility (including its 
very existence) from the irreducibly reversible dynamical 
principle of Standard Hamiltonian mechanics, we must 
complement it with some additional postulate that seems 
to contradict it. 

Within our nonlinear quantum (thcrmo)dynamics 
bascd on Eqs. (|ï|) and (||) the paradox is resolved. Thc 
augmented state domain ansatz broadens the set of con- 
ceivable states but includes the Standard pure states, and 
the nonlinear equation of motion describes irreversibile 
time evolutions and entails reciprocity and fluctuation- 
dissipation relations, but maintains the Standard unitary 
dynamics of pure states. 

Finally, from the Heisenberg time-energy uncertainty 



13 



principle I derived a lower bound for the internal- 
relaxation-time functionals from which follows an up- 
per bound for the rate of entropy production. Con- 
sequently, I proposed a physically intriguing maximal- 
cntropy-production-rate ansatz, by which each indivisi- 
ble subsystem follows the direction of steepest perceived 
entropy ascent at the highest rate compatible with the 
time-energy uncertainty principle. If this ansatz is ex- 
perimentally verified, the nonlinear dynamics is complete 
and self-consistent, with no need of new physical con- 
stants. 
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APPENDIX A: CRITERIA FOR A GENERAL 

(NONLINEAR) QUANTUM DYNAMICS 
COMPATIBLE WITH THERMODYNAMICS 

Within a quantum theory that accepts the augmented 
set of true quantum states described by state operators 
p without the restriction p 2 — p, and a nonlinear dynam- 
ical law for a strictly isolated system, the following de- 
manding set of conditions should be satisfied in order for 
the theory to be compatible or, better, imply the second 
law of thermodynamics without contradicting the fun- 
damental results of Standard quantum mechanics (QM). 
Obviuosly thcse are the criteria I followed in designing 
Eqs. (|ï|) and (g), and discussed at length in Refs. H|). 

Conditions 6, 7 and 8 are closely related to the con- 
dition recently referred to as strong separability (TH) . I 
also added a condition on correlation and entanglement 
to reflect the need to avoid, and possibly resolve, physi- 
cal inconsistències related to nonlocality issues, as well as 
a strong causality condition that is nontrivial and quite 
demanding both from the conceptual and the technical 
mathcmatical points of view. 

Certainly, when viewed from different perspectives — 
e.g., different physical interpretations of the augmented 
state domain p 2 ^ p ansatz, of the role of the nonlinear 
extension of the Schròdinger equation of motion, of the 
Shannon-von Neumann entropy functional — fc B Tr(plnp) 
versus other nonextensive functionals, of the role of the 
system's environment and the measuring apparati, and 
so on — some authors might view this set of conditions 
as too strong in many respeets. Nevertheless, my equa- 
tion of motion demonstrates that at least a satisfactory 
dynamics exists which satisfies all such conditions and, 
in my view, features a number of intriguing, unifying and 
far-reaching implications. 



1. Causality. Forward and backward in time 

Considering the set V of all linear, hermitian, 
nonnegative-defmite, unit-trace operators p on Tí, every 
solution of the equation of motion, i.e., every trajectory 
u(t, p) which at time t = passes through state p in V, 
should lie entirely in V for all times í, — oo < t < +oo. 

2. Standard QM unitary evolution of p 2 — p states 

The unitary time evolution of the states of QM ac- 
cording to the Schròdinger equation of motion must be 
compatible with the more general dynamical law. These 
trajectòries, passing through any state p such that p 2 = p 
and entirely contained in the state domain of Quantum 
Mechanics, must be solutions also of the more general 
dynamical law. In view of the fact that the states of QM 
are the extreme points of the augmented state domain, 
the trajectòries of QM must be boundary solutions (limit 
cycles) of the dynamical law. 

If the complete dynamics preserves the feature of 
uniqueness of solutions throughout the augmented state 
domain, then pure states can only evolve according to 
the Schròdinger equation of motion. In this case, no tra- 
jectory can enter or leave the state domain of QM and 
by continuity, therc must be trajectòries that approach 
indefmitely these boundary solutions (of course, this can 
only happen backward in time, as t — > — oo). 

3. Conservation of energy and number of partides 

If the system is isolated, the value of the energy func- 
tional e(p) — Tr(pií), where H is the Standard QM 
Hamiltonian operator, must remain invariant along ev- 
ery trajectory. If the isolated system consists of a vari- 
able amount of a single type of particle with a number 
operator N that commutes with H, then also the value 
of the number-of-particle functional n(p) = Tv(pN) must 
remain invariant along every trajectory. If the isolated 
system consists of r types of partides each with vari- 
able amount and each with a number operator N that 
commutes with the Hamiltonian H , then also the value of 
each number-of-particle functional rii(p) = Tr (pNi) must 
remain invariant along every trajectory. Depending on 
the type of system, there may be other time-invariant 
functionals. 



4. Stability of the thermodynamic equilibrium 
states. Second law 

A state operator p represents an equilibrium state if 
dp/dt = when the system is isolated [e.g., H ^ H(t)]. 
For each given set of feasible vàlues of the energy func- 
tional e(p) and the number-of-particle functionals rii(p) 
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(i.e., the functionals that must remain invariant accord- 
ing to Condition 2 above), among all the equilibrium 
states that the dynamical law may admit there must 
be one and only one which is globally stable (definition 
and discussion in Appendix [b]) . This stable equilibrium 
state must represent the corresponding state of equilib- 
rium thcrmodynamics and, therefore, must be of the form 
given by Eq. (£9j). All the other equilibrium states that 
the dynamical law may admit must not be globally sta- 
ble. 



5. Entropy nondecrease. Irreversibility 

The principle of nondecrease of entropy must be sat- 
isfied, i.e., the rate of change of the entropy functional 
— fc B Tr(plnp) must be nonnegative along every trajec- 
tory. 



6. Non-interacting subsystems. Separate energy 
conservation 

For an isolated system composed of two distinguish- 
able subsystems A and B with associated Hilbert spaces 
TÍ A and 7í B , so that the Hilbert space of the system is 
H = H A ®Tt B , if the two subsystems are non-interacting, 
i.e., the Hamiltonian operator H = Ha<E)Ib + Ia®Hb : 
then the functionals Tt[(Ha i 8>Ib)p] = ^a(HaPa) and 
Tï[(Ia®Hb)p] = Tïb(HbPb) represent the energies of 
the two subsystems and must remain invariant along ev- 
ery trajectory, even if the states of A and B are corre- 
lated, i.e., even if p ^ pa®Pb- Of course, pa — Trs(p), 
Pb = Tta(p), Tr^ denotes the partial trace over 7í B and 
Tr^ the partial trace over H. A . 



7. Independent states. Weak separability. Separate 
entropy nondecrease 

Two distinguishable subsystems A and B are in inde- 
pendent states if the state operator p = pa®Pb, so that 
the entropy operator S = — fc B lnp = Sa^Ib + -Oi ( S><S',b = 
—k B [BA ln pa®Ib+Ia®Bb hips]. For permanently non- 
interacting subsystems, every trajectory passing through 
a state in which the subsystems are in independent states 
must maintain the subsystems in independent states 
along the entire trajectory. When two uncorrelated Sys- 
tems do not interact with each other, each must evolve 
in time independently of the other. 

In addition, if at some instant of time two subsystems 
A and B, not necessarily non-interacting, are in inde- 
pendent states, then the instantaneous rates of change of 
the subsystem's entropy functionals — k B Tv(pA hi pa) and 
— fc B Tr(pB h\ps) must both be nondecreasing in time. 



8. Correlations, entanglement and locality. Strong 
separability 

Two non-interacting subsystems A and B initially in 
correlated states (possibly due to a previous interaction 
that has then been turned off) should each proceed in 
time towards less correlated states or, at least, maintain 
the same level of quantum entanglement. The genera- 
tion of quantum entanglement between interacting sub- 
systems should emerge only through the Schròdinger-von 
Neumann term —i[H, p]/h of the equation of motion, 
whereas the other terms, that might entail loss of cor- 
relations between subsystems, must not be able to create 
them. This condition is perhaps too strong and its valid- 
ity for my equation is still only conjectural. In any case, 
the dynamics should not generate locality problems, i.e., 
faster-than-light communication between noninteracting 
subsystems, even if in entangled or correlated states. In 
other words, when subsystem A is not interacting with 
subsystem B, it should never be possible to influence the 
local observables of A by acting only on the interactions 
within B, such as switching on and off parameters or 
measurement devices within B. This does not mean that 
cxisting correlations between A and B established by past 
interactions should have no influence whatsoever on the 
time evolution of the local observables of either A or B. 
In particular, I see no physical reason to request that 
two different states p and p' such that p' A = pa should 
evolve in such a way that àp' A /àt = dp^/dí whenever 
A is isolated (but not uncorrelated) from the rest of the 
overall system. For example, state p' could be the max- 
imum entropy stable equilibrium state (and, therefore, 
p' = p'a®Pb' dp'/dí = 0) whereas in state p subsystems 
A and B could be correlated and evolving in time to- 
wards the stable equilibrium state or p could even be a 
pure entangled state evolving along a unitary trajectory 
according to Condition 2 above, and therefore it would 
never reach stable equilibrium. 

APPENDIX B: LYAPUNOV STABILITY AND 
THERMODYNAMIC STABILITY 

The condition concerning the stability of the thermo- 
dynamic equilibrium states is extremely restrictive and 
requircs furthcr discussion. 

In ordcr to implement Condition 5 in Appendix |A|, 
we need to establish the relation between the notion of 
stability implied by the second law of thermodynamics 
pp| , pï| and the mathematical concept of stability. An 
equilibrium state is stable, in the sense required by the 
second law, if it can be altered to a different state only 
by interactions that leave net effects in the state of the 
enviromment. We call this notion of stability global sta- 
bility. The notion of stability according to Lyapunov is 
called local stability. 
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We denote the trajectòries generated by the dynamical 
law on our state domain by u(t,p), i.e., u(t,p) denotes 
the state at time t along the trajectory that at time t = 
passes through state p. A state p e is an equilibrium state 
if and only if u(t, p e ) = p e for all times t. An equilibrium 
state p e is locally stable (according to Lyapunov) if and 
only if for every e > there is a S(e) > such that 
d{p, p e ) < 5(e) implies d(u(t, p), p e ) < e for all t > and 
every p, i.e., such that every trajectory that passes within 
the distance 8{e) from state p e proceeds in time without 
ever exceeding the distance e from p e . Conversely, an 
equilibrium state p e is unstable if and only if it is not 
locally stable, i.e., there is an e > such that for every 
6 > there is a trajectory passing within distance S 
from p e and reaching at some later time farther than the 
distance e from p e . 

The Lyapunov concept of instability of equilibrium is 
clearly equivalent to that of instability stated in ther- 
modynamics according to which an equilibrium state is 
unstable if, upon experiencing a minute and short lived 
influence by some system in the environment (i.e., just 
enough to take it from state p e to a neighbouring state 
at infinitesimal distance S), proceeds from then on spon- 
taneously to a sequence of entirely different states (i.e., 
farther than some finite distance e). 

It follows that the concept of stability in thcrmody- 
namics implies that of Lyapunov local stability. How- 
ever, it is stronger because it also excludes the concept of 
metastability. Namely, the states of equilibrium thermo- 
dynamics are global stable equilibrium states in the sense 
that not only they are locally stable but they cannot be 
altered to entirely different states even by means of in- 
teractions which leave temporary but finite effects in the 
environment. Mathematically, the concept of metasta- 
bility can be defined as follows. An equilibrium state p e 
is metastable if and only if it is locally stable but there 
is an 77 > and an e > such that for every 5 > there 
is a trajectory u(t, p) passing at t = between distance 
77 and i] + 5 from p e , T) < d(u(0, p), p e ) < 77 + S, and 
reaching at some later time t > a distance farther than 
77 + e, d(u(t, p), p e ) >?? + £■ Thus, the concept of global 
stability implied by the second law is as follows. An equi- 
librium state p e is globally stable if for every 77 > and 
every e > there is a S(e, 77) > such that every tra- 
jectory u(t, p) with 7] < d(u(0, p), Pe) < r) + S(e, 77), i.e., 
passing at time t = between distance 77 and 77 + S from 
p e , remains within d(u(t, p), p e ) < 77 + e for every t > 0, 
i.e., proceeds in time without ever exceeding the distance 
r] + e. 

The second law requires that for each set of vàlues of 
the invariants Tr(pH) and Tr(pGi) (as many as required 
by the structure of the system), and of the parameters 
embedded in the Hilbert space Tí and the Hamiltonian 
H describing the external forces (such as the size of a 
container) , there is one and only one globally stable equi- 
librium state. Thus, the dynamical law may admit many 
equilibrium states that all share the same vàlues of the 
invariants and the parameters, but among all these only 



one is globally stable, i.e., all the other equilibrium states 
are either unstable or metastable. 

Interestingly, we may use this condition to show that 
a unitary (Hamiltonian) dynamical law would be incon- 
sistent with the second-law stability requirement. A uni- 
tary dynamical law in the augmented kinematics would 
be expressed by an equation of motion ifip = [H, p] 
with trajectòries u(t,p) = E/(t)jO?/ -1 (t) with U(t) = 
exp(—itH/K). Such a dynamical law would admit as 
equilibrium states all the states p e such that p e H = Hp e . 
Of these states there are more than just one for each set 
of vàlues of the invariants. With respect to the mètric 
d{piiP2) — Tr|pi — P2I, it is easy to show that every tra- 
jectory u(t, p) would be equidistant from any given equi- 
librium state p e , i.e., d(u(t, p), p e ) — d(u(0, p), p e ) for all 
t and all p. Therefore, all the equilibrium states would 
be globally stable and there would be more than just one 
for each set of vàlues of the invariants, thus violating the 
second-law requirement. 

The entropy functional — fc B Tr(plnp) plays a uscful 
role in proving the stability of the states of equilib- 
rium thermodynamics [Eq. (|29|)] provided that the dy- 
namical law guarantees that — fc B Tr[?j(í, p) lnit(í, p)] > 
—k B TT(phip) for every trajectory, i.e., provided Condi- 
tion 6 above is satisfied. The proof of this is nontrivial 
and is given in Ref. f2Ql where, however, we also show 
that the entropy functional is not a Lyapunov function, 
even if, in a strict sense that depends on the continuity 
and the conditional stability of the states of equilibrium 
thermodynamics, it does provide a criterion for the sta- 
bility of these states. Anyway, even if the entropy were a 
Lyapunov function, this would suffice only to guarantee 
the local stability of the states of equilibrium thermody- 
namics but not to guarantee, as required by the second 
law, the instability or metastability of all the other equi- 
librium states. 



APPENDIX C: CHARACTERISTIC TIMES 

Using Eq. @ , the rate of change of the mean functional 
f(p) = Tr(pF) may be written as 



d/Q>) 
df 



2(y/pF\E) 



(Cl) 



For the Schròdinger-von Neumann evolution, the char- 
acteristic time of change of f(p) may be defined as (^?J 



1 _ [d/( P )/dí] 2 _4(^pAF\E H ] 



'FH 



(AFAF) 



(AFAF) 



(C2) 



Because operators ^fpAFj^J (AFAF) are unit norm, 
in the sense that {^fpAF\jpAF)/(AFAF) = 1, it fol- 
lows that the characteristic times tfh are bounded by 
the value attained for an operator F such that ^/pAF is 
in the same direction as E H , i.e., such that 
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E 



H 



Therefore, 



<±{ E h\ E h, 



(C3) 



(C4) 



FH 



' H 



For this reason, in Eq. ( |Ï4| ) we take th equal to the lower 
bound of the tfh's- 

By analogy, but considcring the understanding of the 
contribution of each subsystem to the overall-system dy- 
namics embedded in our nonlinear dynamics, we define 
the characteristic time of the dissipative change of the 
mean functional f(p) = Tr (fiF) due to the J-th con- 
stituent subsystem as 



1 



[d/( P )/dí] 2 _ 4 (vpj( AF y \ e dj 



FD.J 



(AFAF) J 



(AFAF)' 7 



(C5) 



For such special choice of the Xj's, assuming X\ = /, 
we have x 1 (p) = Tr(p) = 1, x^i(p) = Tr(pXj^i) = 0, 



(AX 1 AX i ) = 0, (AX m AX m ) 



Lij(p) = and 



Dxi(p)/Dt = for i or j < a, Lij(p) = 5ij/r(p) and 
Dxi(p)/Dt = fi(p)/r(p) for i and j > a, so that 



(D2) 
(D3) 
(D4) 




dí 



y/(AHAH) J£f k {py< 



k>a 



where in the last two equations we made use of Eqs. ( |25| ) 
and Relations ( |53| ) and (|^) with strict equality. 

In view of the dependence of the Xj 's on p and, t hcrc - 
fore, on time, the apparently simple form of Eq. ( pl| ) 
may not be as useful as it seems. 



Again, because operators v /p~j(AF) J / y (AFAF) are 
unit norm, in the sense that 

(^fp- J {AF) J \^fp- J {AF) J )/(AFAF) J = 1, it follows that 
the characteristic times t fdj are bounded by the value 

attained for an operator F such that ^/]TJ(AF) J is in the 
same direction as E DJ , i.e., such that 



/P7(AF) J 
(AFAF)' 1 



E 



D.J 



E D,J \ E DJ, 



Therefore, 



<H E Dj\ E DJ, 



FD.J 



1 

DJ 



(C6) 



(C7) 



For this reason, in Eqs. (p0|) and (10S) we take tjj and 



t dj ec l ua l to the respective lower bounds of the r 



FD.J 



APPENDIX D: SPECIAL FORM OF THE 
EQUATION OF MOTION 

If, for a given state operator p, we construct the set 
{y/pXj} so as to be an orthogonal extension of the 
orthonormal subset {^/pAi}, i.e., with Xj = Ai for 
i < a and (^fpX^^fpXj) — for all i and j, then 
D = — J2j> a fj y/P-X-j [Eq. (|3^)] and the onc-constituent 
equation of motion reduces to the (only apparently lin- 
ear) form 



dp 
dí 



7; 



[tf.d + aT^E/iÍpH^OO.p}. ( D1 ) 

T ^ P > j>a 



where the dependences of fj , Xj and t on p are evidenced 
in order to emphasize the nonlinearity. 



APPENDIX E: EXTENSION TO TIME-VARYING 
RATES OF THE GENERATORS OF THE 
MOTION 

For the purposes of quantum thermodynamics, in my 
view the dissipative part of the equation of motion should 
not account for rates of change of the mean vàlues of the 
generators of the motion other than through the Hamil- 
tonian term, consistently with all the results of Standard 
QM. 

Nevertheless, the mathematical extension of Eqs. @) 
and (||) to (artificially) imposed rates fj ^ (arbitrarily 
specified as functions of time) is straightforward and may 
be useful in applications or other frameworks H. 

It does emerge naturally from the maximization prob- 
lem (p6|); it does so, however, implicitly, through sub- 
stitution of Eq. ( |58| ) back into the constraint Equations 
(fÏ6|b). Instead, the explicit form in terms of projections 
and, therefore, the equivalent expressions by means of 
Gram determinants amount to assuming, for the opera- 
tor D in the equation of motion, instead of Eq. (j^c) or 
the equivalent Eqs. (|2C|), 



D = yp~\a.p]±c{^pRi} 



(El) 



where the explicit expression of the j-th term in the sum- 
mation is given by Eq. (po|b) where we interchange ev- 
ery where y/p ln p with y/pRj and 



T[P,r] T{Jp]n P ,{JpIU}) 



(E2) 



Geometrically, the additional terms are in the steepest- 
rj-ascent direction compatible with maintaining constant 
the mean vàlues of the other generators and the entropy. 
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Analogous obvious extension to the composite system 
case can be obtained by assuming for the operator Dj, 
instead of Eq. ([[]) , 



formally better, and suggest the apparently alternative 
equation of motion based on the following definitions 



D.J = [\/P7(^lnp) J ] ±£{v77 - (fli) , 



n > (E3) 



with the ctj's such that 



APPENDIX F: EXTENSION TO OTHER 
ENTROPY OR MEAN VALUE FUNCTIONALS 

For the purposes of quantum thermodynamics, in my 
view the necessary entropy hmctional is — fc B Tr(plnp), 
for the reasons in Refs. . 

Howcver, in view of the recent literature on nonextcn- 
sive quantum theories, as suggested also in Ref. ]Ï2| ], it 
may be useful to note that the entire formalism of my 
equation of motion can be readily reformulated in the 
case of any other well-behaved entropy functional s(p) 
and set of nonlinear generator functionals (p) . 

For the single constituent system it sumccs to sub- 
stitute throughout the operator (— 2k By /pI+)y/pS — 
(—2k B ^/~p) — 2fc B y / plnp (notice that in the Gram deter- 
minants the addenda in paranthescs cancel out) with the 
new entropy gradient operator, ds(p)/dy/p, and the op- 
erators 2^/pRi with [operators which when symmetrized 
({A, A^}/2) are equal to] the gradient operators of the 
new generator functionals, dri(p) / d^fp. 

However, for a composite system, consistently with the 
nonextensivity of these theories, it may be difficult or not 
at all possible to identify the operators corresponding to 
(S) J and (Ri) J representing the subsystems' local per- 
ceptions of the entropy and the generators. 



APPENDIX G: A NOTEWORTHY EQUATION 
FOR A COMPOSITE SYSTEM THAT FAILS TO 
MEET A SEPARABILITY CONDITION 

It is interesting to note that the role of the ^fp~j oper- 
ators in Eq. (0) [and that of Jp~ in Eq. (|J)] is formally 
useful but only auxiliary, because wherever there is a 
yfpj, another ^/p~J comes in front or behind it. In fact, 
in my doctoral thesis Q the equation of motion was in 
terms of Pj only. I realized the usefulness of the ^/pj 
formalism only later, in connection with the proof of the 
steepest-entropy-ascent geomètric property |J . 

By allowing a more relevant role of the yíp operator, 
an alternative to the construction in Eq. may appear 



E T 



L·j=i 



i ry 



{^-pF)'l = Tr 7 [{I^ 



f>7> 



D[ r = [(^p\n P y v 



(Gla) 
(Glb) 
(Glc) 



In fact, it can be readily verified that for a given 
F on Tí with [F, H] = the rate of change of 

Tr(pF) would be zero if and only if (y/pí 1 )^ is in 

£{(VP)^ (VpH) j v , [(y/pdyj}, and so Tr(p), Tr(pH), 
and Tr(pGí) would be conserved. The expressions for the 
rate of entropy production, the Onsager relations and the 
other results would be almost identical to those obtained 
from my equation, except for the substitution throughout 

of ^(íy with (^pFy. 

However, the resulting dynamics would fail to satisfy 
at least the important property expressed by Eq. (|6^), 
because it can be verified that if H = H + I j<E>Hj 
but p pj®p-j then separate conservation of the non- 
interacting subsystem's energy would not be guaranteed. 
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